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ABSTRACT 



The delayed detonation model describes the observational properties of the 
majority of type la supernovae very well. Using numerical data from a three- 
dimensional deflagration model for type la supernovae, the intermittency of the 
turbulent velocity field and its implications on the probability of a deflagration- 
to-detonation (DDT) transition are investigated. From structure functions of the 
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turbulent velocity fluctuations, we determine intermittency parameters based on 
the log-normal and the log-Poisson models. The bulk of turbulence in the ash 
regions appears to be less intermittent than predicted by the standard log-normal 
model and the She-Leveque model. On the other hand, the analysis of the tur- 
bulent velocity fluctuations in the vicinity of the flame front by Ropke suggests a 
much higher probability of large velocity fluctuations on the grid scale in compar- 
ison to the log- normal intermittency model. Following Pan et al., we computed 
probability density functions for a DDT for the different distributions. The de- 
termination of the total number of regions at the flame surface, in which DDTs 
can be triggered, enables us to estimate the total number of events. Assuming 
that a DDT can occur in the stirred flame regime, as proposed by Woosley et 
al., the log-normal model would imply a delayed detonation between 0.7 and 0.8 
seconds after the beginning of the deflagration phase for the multi-spot ignition 
scenario used in the simulation. However, the probability drops to virtually zero 
if a DDT is further constrained by the requirement that the turbulent velocity 
fluctuations reach about 500 km s -1 . Under this condition, delayed detonations 
are only possible if the distribution of the velocity fluctuations is not log-normal. 
From our calculations follows that the distribution obtained by Ropke allow for 
multiple DDTs around 0.8 seconds after ignition at a transition density close to 
1 x 10 7 gcm~ 3 . 

Subject headings: stars: supernovae: individual: la — hydrodynamics — turbu- 
lence — methods: statistical 



Introduction 



Although deflagration mode 



s of Type la supernovae (e.g. iReinecke et al.ll2002l ; iGamezo et al 



199ll ; iRopke fc Hillebrandtl 120051) reproduce bulk prope rties of fainter Type la supernovae 



delayed detonations (see iHillebrandt &: Niemeyerl |2000| . for a review of explosion scenarios) 
seem to be the the most promising way of m odeling the majority of the observed events 
faopke & Niemeverll2005l : iMazzali et alll2007h . 



However, a theoretical explanation for the deflagration-to-detonation (DDT) transition 
that gives rise to this type of ther monuclear explosion has be en missing so far. Originally 
postulated by iKhokhlovl (119911 ) and IWooslev &: Weaverl (119941 ). the possibility of a DDT in 
type la supernovae was questioned by lNiemeyer Jl999h . because the preconditioning required 
to trigger a self-sustaining detonation front would be extremely unlikely to occur i n the defla- 



gration phase. This conclusion was based on ensemble-average scaling arguments. iPan et al. 
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( 120081 ). on the other hand, pointed out that highly intermittent turbulent fluctuations act in 
favour of a DDT, because intermittency allows for velocities much higher than the ensemble 
average. 



As proposed by lNiemeyer fc Woosleyl (119971 ). iNiemeyer fc Kersteinl ( 119971 ) and lKhokhlov et al 



(119971 ). a necessary condition for trigge ring a DDT is tha t the distributed burning regime 
has t o be entered. Based on log-normal (jKolmogorovlll962l ) and log-Poisson (jShe fc Leveque 
19941 ) intermittency models with typical parameters expected for turbulence in thermonu- 
clear supernovae, Pan et al. calculated probabilities for this DDT condition. Their calcula- 
tions i mplied transition dens ities in the ran g e 3.8 x 10 7 . . . 2.7 x 10 8 gcm -3 . iLisewski et al. 
( 120001 ) . IWoosleyl (120071 ). and IWoosley et al.l (120091 ) , on the other hand, provided numerical 
evidence that there are stronger constraints on a DDT than just entering the regime of 
distributed burning. As a consequence, lower transition densities would be favored. Analyz- 
ing the small-scale turbulent velocity distrib ution p r edicte d by the subgrid-scale model in 
several thermonuclear supernova simulations, iRopkd (120071 ) found that intermittent velocity 
fluctuations up to 1000 km s -1 on a length scale of 10 km are possible. Fluctu ations of this 



magnitude could trigger a detonation in accordance with ILisewski et al.l (120001 ) 



The ambiguity of conditions for DDTs in type la supernovae makes predictions difficult. 
Here, we follow a pragmatic approach and investigate different DDT criteria, assuming sim- 
ple approximations to the microphysical conditions. Following the approach of Pan et al., we 
calculate the proba bility of a DDT usi n g data sets from a highly resolved three-dimensional 
deflagration model (IRopke et al.l 120071 ) . ICiaraldi-Schoolmann et al.l (120081 ). in the following 
referred to as paper I, determined characteristic scales and scaling exponents from the struc- 
ture functions of the turbulent velocity fluctuations up to order six. We fitted log-normal 
and log-Poisson intermittency models to these scaling exponents. In addition, we determined 
probability density functions (pdfs) of the mass density constrained to the flame surface in 
order to calculate the probability of a DDT at different stages of the deflagration. There are 
no further adjustable parameters apart from the initial conditions of the deflagration model. 
Assuming a certain size for DDT regions, the total number of DDTs can be obtained from 
the cumulative probability and the number of regions in the volume occupied by the flamed 



Basically, iPan et al.l (120081 ) assumed that a DDT is triggered once the Karlovitz number 
Ka (the square root of the ratio of t he laminar flame width to the Gibson scale as d efined 
in Section [3]) exceeds a certain value (INiemeyer fc Kersteinl 1 19971 ; I Woosley et al.ll2009l ). i. e., 



1 Down to the Gibson scale, the flame front can be considered as a fractal surface, which is numcerically 
represented by a zero level set. Thus, the volume occupied by the flame is defined by all grid cells within a 
certain distance to the zero level set. 
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burning has advanced sufficiently far into the regime of distributed burning. Because the 
critical size of a preconditioned region is extremely sensitive to the mass density, it is difficult 
to estimate the total number of DDTs. Nevertheless, our results suggest that a delayed 
detonation would result at an early stage of the explosion if the conditions for a DDT were 
constrained by the Karlovitz number only. A further implicatio n would be a tra nsition 
density well above 1 x 10 7 gcm~ 3 , which was found in the study by lPan et al.l (120081 ). 



Whereas the criterion f or the onset of distrib uted burning can be expressed in terms 
of laminar flame properties, IWoosley et al.l (120091 ) argued that the conditions for a DDT 
cannot be parameterized on the basis of the laminar flame speed and th ickness. Using 



numerical data for the nuclear time scale in the so-called stirred flame regime (jKersteinll200ll ) 
and following the same procedure as outlined above, we calculated the DDT probability. 
Although there are substantial uncertainties in this approach, a DDT at later time and lower 
density compared to the estimates following from laminar flame properties is implied. Most 
important, we find that a DDT can be excluded if the magnitude of the turbulent velocity 
fluctuati ons is required to exceed about 500kms _1 , which is suggested by microphysical 
studies ( jWoosleyl 120071 ) . Delayed detonations are possible, however, if the distribu tion of 
the ve locity fluctations significantly deviates from a log-normal shape. According to iRopke 
( 120071 ). this appears to be the case for regions near the flame front. 



The paper is organized as follows. In the next Section, we briefly review the log- 
normal and log-Poisson intermittency models. Moreover, the results on the intermittency of 
turbulence in the simulation are discussed. The procedure of calculating the DDT probability 
is described in Section |3j In Section HJ the estimation of the total number of DDTs is 
explained. We present our numerical results in this Section [5j followed by a discussion in 
the last Section. 



2. Intermittency 



According to the Kolmogorov theory of statistically stationary, homogeneous, and isotropic 
turbulence, the structure functions S p (£) := (5v p (£)), i. e., the ensemble average of the veloc- 
ity fluctuations to the power p over a length scale £, follow power laws in the inertial subrange. 
The scaling exponents are given by ( p = p/3. In several experiments, however, departures 
from the p/3 scaling were found for higher-or der structure functions. This phenomenon is 
attributed to the intermittency of turbulence (jFrisch 1995 . see). There are various theoret- 
ical intermittency models that predict anomalous scaling exponen ts. One example is the 
Kolmogorov- Oboukhov model (jKolmogorovl Il962l ; lOboukhovl Il962l ). which assumes a log- 
normal probability distribution function for the turbulent energy dissipation q on a length 



- 5 - 



scale 



pdf(ln( Q /e)) 



I'KO", 



exp 



[ln( Q /e) + af /2]' 
2rf 



(1) 



where e is the mean rate of energy dissipation and cr| = fj,ln(L/£) for the integral scale L. 
The scaling exponents resulting from this model are given by 



/■On) 



|-^(p-3) 



(2) 



where the parameter p, is defined by p = 2 — £ 6 . Experiments and simulations indicate 
^6 ~ 1.8 and therefore p ~ 0.2. A defficiency of the log- normal model is that it fails to 
describe the observed higher-order scalings for p > 10 (see iPan et al.l (120081 ) for a detailed 
discussion of the limi tations of this model). In this respect, the intermittency model by 
She fc Levequd (119941 ) has been particularly successful: 



a(SL) _ P , 2 



1 



p/3' 



(3) 



Dubrulld (119941 ) and IShe fc Waymird (119951 ) showed that this model can be accommodated 
within a two-parameter family of intermittency models, which are based on log-Poisson 
statistics: 



(IP) 



:i -A)f 



A 



1-P 



(4) 



The random cascade factor (3 G [0, 1] specifies the degree of intermittency. Non-intermittent 
Kolmogorov scaling is obtained in the limit (5 — > 1. The parameter A specifies the scaling 
properties of the most intermittent dissipative structures and is related to the codimension 
C of these structures: C = A/(l — 0). Setting A = 2/3 and (3 = 2/3, the She-Leveque 
model ([3]) follows from equation (j4]). For these parameters, the codimension is C = 2. 
The interpr etation is that vortex filaments are the most dissipative structures in subsonic 
turbulence (IShe fc Levequd I1994T) . For superson i c turbulence, 1 < C < 2 because shock 



fronts also dissipate energy (IBoldyrev et al.ll2002l ; iKritsuk et al.l 120071 : ISchmidt et al.ll2009l ). 



Equation (j2J) as well as imply ^ 3 = 1 and, hence, Ss(£) ~ e£. In paper I, we 
determined the scaling exponents form power-law fits to the structure functions for p < 6. 
It was found that there is a transition length £k/rt ~ 10 km, where the radial second- 
order structure function changes from Kolmogorov scaling (S2(£) oc £ 2 ^ 3 ) to Rayleigh- Taylor 
scaling (S^t) oc £). The values of £k/rt for t > 0.5 seconds are listed in Tabled! While £3 
is approximately unity for the angular structure functions and also for the radial structure 
functions in the subrange £ < £k/rt, the radial third-order exponent for £ > £k/rt is about 
1.5. On the other hand, it was shown in paper I that the relative scaling exponents Z p = 
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Cp/C3 of the radial structure functions are nearly equal for the two subranges. However, 
since a consistent treatment of an exponent £3 different fro m unity is highly non-trivial in 



the fra mework of the intermittency models mentioned above (jDubrullelll994l ; ISchmidt et al. 



2008bl ). we do not consider the radial scaling exponents for length scales greater than £k/rt- 



The values of /x resulting from fits of the log-normal model to the exponents ( p for t = 
0.5, ... ,0.9 seconds are summarized in Table [TJ Table [2] demonstrates that these values are 
consistent with // = 2 — £ 6 within the error bars of the scaling exponents. For t = 0.7 seconds, 
the scaling exponents of the radial and angular structure functions are plotted together with 
the fit functions in Figure [TJ For comparison, also the log-normal model with fi = 0.2 and 
the She-Leveque model ([3]) are plotted. We find \x ~ 0.1. Clearly, the numerical data from 
the supernova simulation indicate a lower degree of intermittency than these two models. 
Moreover, there is a remarkable agreement between the radial scaling exponents for the range 
of length scales t < £k/rt and the scaling exponents of the angular structure functions over 
the whole range of length scales (also see paper I). 

Fitt ing the ge n eral l og-Poisson model ([3]) to the scaling exponents, yields ambiguous 
results. IPan et al.l (120081 ) propose to set A = 2/3. This appears to be a sensible choice 
on grounds of the hypothesis of a universal dissipation time scale, which determines A. 
However, the best fits constrained by A = 2/3 imply a codimension greater than 2. Since 
it is hard to come up with a reasonable physical interpretation of this result, we tested the 
hypothesis C = 2, i. e., we consider the codimension of the most intense dissipative structure 
to be fixed. The resulting fit functions are almost indistinguishable from the best log-normal 
fits (see Figured]). The log-Poisson fit parameters are also listed in Table [TJ As one can 
see, (3 > 2/3, indicating a lower degree of intermittency in comparison to the She-Leveque 
model, and A m 0.5. In this respect, the trend indicated by the log-normal model fits is 
confirmed. The anomalous intermittency parameters point at deviations from fully developed 
turbulence. A possible reason is that turbulence in supernova explosions does not reach a 
statistically stationary st ate. Since there are indicati ons for A > 2/3 in the case of highly 
compressible turbulence (ISchmidt et al.ll2008bl . 120091 ). the lower value of A in the present 
case cannot be attributed to compressibility. Since the log-Poisson model is sensitive even to 
small errors in the higher-order scalings, the deviations from the She-Leveque model might 
be spurious. In this regard, however, one has to keep in mind that insufficiently converged 
statistics us ually causes an overestimate of intermittency rather than an underestimate (see 



Frischlll995l ). 



Assuming that the low value of /i is a peculiar property of non-stationary RT-driven tur- 
bulence in thermonuclear supernovae, one might raise the question if the small-scale isotropy 
discussed in paper I is genuine or, possibly, an artifact of the employed SG S model. Since 



our SGS model features a localized closure for the turbulence energy flux (ISchmidt et al. 
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20061 ). isotropy on numerically resolved scales is not required. Only the unresolved turbulent 
velocity fluctuations are assumed to be isotropic. Fig. 2 in paper I demonstrates that the 
resolved velocity field is isotropic on length scales smaller than £k/rt- At late time, when 
£k/rt becomes less than the size of the grid cells, the radial and angular structure functions 
deviate even on the smallest numerically resolved scales and, thus, we have anisotropy on 
these scales. Consequently, isotropy does not appear to be forced by the SGS model if there 
is a pronounced anisotropy on length scales near the numerical cutoff scale. Although we 
cannot exclude the possibility that the SGS tends to make the nearby resolved scales more 
isotropic than in reality, because the interactions between resolved and unresolved scales are 
not strictly local, it is plausible that the observed isotropy on length scales £ < £k/rt is 
basically genuine. 

Analyzing the statistics of the SGS velocity fluctuations in the vicinity of the flame 



front. iRopkd (120071 ) found the high- velocity tail of the pdf is very well fitted by the expression 
pdf(8v(£)) = exp [ao5v ai (£) + 02] f° r Sv(£) > 10 7 cms _1 , where £ = A. The local velocity 
fluctuation at the numerical cutoff is given by the specific SGS turbulence energy: Sv(A) = 



Qsgs = A/2A; S g S (jSchmidt et al.ll2006l ). The coefficients a , oi, and a 2 are fitted for each instant 
separately. The variation of the coefficients is due to the evolution of the turbulent flow. 
Changes might also result from the temporal shift of the cutoff scale. For the phase of the 
explosion we are interested in, A is about 10 km, which in turn is close to £k/rt- Further 
studies indicated that the above form of the pdf also holds for resolved velocity fluctuations on 
length scales £ ~ £k/rt- Thus, we assume that the fit coefficients do not change significantly 
for nearby scales at a fixed time. The corresponding pdf for q = 5v 3 (£)/£ is given by 

£1/3 

pdf (e e ) = — ^ exp [a (£e e ) ai/3 + a 2 ] for £ ~ 10 km. (5) 



Fig. 1 in IRopkd (120071 ) shows that the high- velocity tail is much flatter than the prediction of 
the log-normal intermittency model. Consequently, there appears to be a much higher degree 
of intermittency. It is important to note that this applies to velocity fluctuations at the flame 
surface. In this regard, the behavior of turbulence close to the flame front appears to be 
markedly different from the properties of the bulk of turbulence in ash regions, for which the 
log-normal fits to the scaling exponents imply a relatively low degree of intermittency. We 
will see in Section [5] that this difference is crucial for the occurrence of delayed detonations 
if tight constraints on DDTs are assumed. 
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3. Deflagration-to-detonation transition probability 



Turbulent deflagration in t he flamelet regime is characterized by a flame thickness £r < 
£q (INiemeyer fc Kersteinl Il997l ). where the Gibson scale £q is the length scale for which 
the turbulent velocity fluctuations are comparable to the laminar flame speed: 5v{£q) ~ 
•5iam- According to the refined similarity hypothesis of iKolmogorovl (119621 ). we have Sv(£) ~ 



e 1 / 3 ^ 1 ' 3 . Assuming isotropic turbulence, the mean rate of dissipation e = V 3 /L, where 
L and V are the integral length and the associated velocity, respectively. However, it was 
shown in paper I that turbulence in thermonuclear supernovae is anisotropic on length scales 
considerably smaller than L. We will therefore give a more precise definition below. In the 
ensemble average, it follows from £q = (s 3 m / V") 3 L that e ~ Ka 2 6fl, where := s 3 am /^fl and 
the Karlovitz number Ka = (^a/^g) 1 ^ 2 ( Petersll2000). Once £g becomes sma l ler than Ia, i. e.. 



Ka > 1 , distributed burning will commence (INiemeyer fc Woosleyl 119971 ; iKhokhlov et al. 



19971 ). It should be noted that the meaning of Ka in terms of the Gibson length becomes 
elusive in the limit Ka 3> 1, because the flame thickness and the laminar burning speed are 
not well defined in this regime. 



3.1. Estimation of the DDT probability based on laminar flame properties 



In log-normal as well as log-Poisson intermittency models, the rate of dissipation av- 
eraged over a region of size I is defined by a random variab le q with a certain probability 
density function pdf(e^). As pointed out by lPan et al.l (120081 ). the probability of distributed 
burning with a certain Karlovitz number in a region of size I is then given b}@ 



/•oo 

P{e t > Ka 2 e fl ) = / pdf (ej) de' e . 

>/Ka 2 e fl 



(6) 



As the Karlovitz number, for which a DDT can occur, increases, the probability given by 
the above equation declines. Apart from Ka, the threshold for eg is given by laminar flame 
properties, because ea is defined by the speed and width of laminar flames. 

Assuming that a DDT is caused by the Zel'dovich mechanism, there must be a sufficient 
number of preconditioned regions that reach the critic al size £ c . This length scale depends 
on the mass den sity and the composition of the fuel (INiemeyer fc Woosleyl 119971 ). While 
Pan et al.l (120081 ) assumed a fixed density and a spherical flame, we evaluate the simulation 
data described in paper I. Specifically, we compute the effective probability that a single 



2 They use K = Ka 2 ^ 3 as a fudge parameter, without referring to the Karlovitz number. 



- 9- 



DDT occurs anywhere near the flame front from the convolution of the conditional proba- 
bility P(ee > Ka 2 efl) with the normalized probability density function of the mass density p 
constrained to the flame front: 



POO 

P DDT (Ka) = / P{e tc > Ka 2 e fl ) pdf (p\G = 0) dp. 
Jo 



(7) 



The flame front is numerically defined by the zero level set, G — 0, where G denotes the the 
level set function. Basically, G(t, r) is a distance function that is positive inside the flame 
(in ash regions) and negative outside. For the numerical calculation of pdf(p|G = 0), the 
grid cells in which G(t,r) switches sign are identified. The distribution of the mass density 
over the flame surface is mostly due to the density stratification in the exploding star. Since 
turbulence in SNe la is weakly compressible, we ignore correlations between local fluctuations 
of the density and the velocity. 

For the calculation of -P(q c > Ka 2 efl), the probability density function of ei has to be 
modelled. From the log-nor mal probability density function (OQ), the following expression for 



P(e tc > Ka 2 e fl ) is obtained flPan et all 120081 ): 



p/ ^ v 2 \ 1 f ( ln(Ka 2 e fl /e) a lc 
P [e lc > Ka e fl ) = -erfc 



2 y y/2a tc 2^2 / 

For the log-Poisson model, on the other hand, it is not possible to derive an analytic expres- 
sion for P (e c > Ka 2 efl), because the Poisson distributions contributing to pdf(e^) are not 
explicitly known. For this reason, we will not consider the log-Poisson intermittency model 
in the following. 



The restrictions of the log- normal model were discussed at length by iPan et al.l (120081 ). 
Data from high-resolution simulations of turbulence suggest that the log-normal model yields 
a good approximation to the distribution of ei within the 5ag wings. However, it is not clear 
whether this range also applies to turbulence in thermonuclear supernova. A conservative 
estimate can be made as follows: Given a random variable X, the largest X- values, for 
which the probability density function pdf (X) is constrained by the moments of order < p, 
is indicated by the peak of X p pdf(X). Since e<? ~ 5v 3 (£)/£, we have 5v p (£)pdf (Sv(£)) ~ 
5v p ^ 1 (£)pdf (\n[Sv 3 (£)/£e\), where the pdf on the right-hand side is defined by equation (CQ). 
We computed the scaling exponents of the structure functions up to sixth order. For £ = 
10 km and the numerical parameters listed in Table [Q Sv 6 (£)pdf (5v (£)) peaks for velocity 
fluctuations around lOOkrns -1 , which corresponds to a range within 2a e from the maximum 
of the log-normal Q-pdf. Consequently, the model extrapolates the central part of the pdf, 
which can be inferred from the known structure functions, to the far tail. It is important 
to note that the range, for which the log-normal pdf is not well constrained by the scaling 
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exponents, is about the rang e of the asymp totic tail (jSJ) resulting from the analysi s of small- 



scale v elocity fluctuations by iRopkd (120071 ). In order to compare to the results of |Pan et al. 



( 120081 ). however, we will also consider a log-normal shape of the pdf up to 5<r^. 



For the determination of the mean dissipation rate e, it is important to bear in mind that 
intermittency models such as the log-normal model apply to statistically isotropic turbulence. 
It was shown in paper I, that Kolmogorov scaling applies on length scales £ < £k/rt- For 
this reason, it is not correct to set e = V 3 /L, where V and L are the scales of energy 
injection, in the case of RT-driven turbulence. In fact, the rate of dissipation is fixed by 
the transition length £k/rt and the velocity scale Sv(£k/rt)- Since the third-order structure 
function Ss(£) ~ e£, we define e ra( j := 5 , 3 !ra d(^K/RT)/^K/RT- However, one can also set e ang := 
S^ang (L)/L, because the scaling of the angular structure function is unique and consistent 
with the Kolmogorov theory for the whole range of length scales. Table [T] gives an overview 
of the numerical values of e ra d, e ang , and L as functions of time. While both definitions yield 
similar results for the mean rate of dissipation, it is obvious from Figure 3 in paper I that 
e <C V 3 /L ~ 53,rad(-£) / L if V is taken to be the characteristic velocity of the Rayleigh- 
Taylor-driven turbulent flow at the length scale L. On the other hand, it is essential to set 
the integral length scale for the distribution of q to L rather than £k/rt, because intermittent 
velocity fluctuations occur on all length scales £ < L. Although turbulence is anisotropic 
for £k/rt ^ £ ;$ L, we assume that the log-normal model can be applied to the whole 
range of length scales on the basis of the dissipation rate e ang ~ e ra d, because the angular 
structure functions continue to follow the scaling laws of isotropic turbulence for £ > £k/rt- 
In Section we will carry out numerical calculations for both the radial and the angular 
parameter sets. 



3.2. Estimation of the DDT probability in the stirred flame regime 



So far, we h ave assumed that a D DT is only constrained by a value of Ka greater than 
unity. However, IWoosley et al.1 (120091 ) argued that an additional requirement for a DDT is 
that Da > 1, where the Damkohler number Da = T/T nur is th e ratio of a dynamical time 
scale T to the nuclear burning time scale r nuc ( Kerstein 200 lh. This regime is called the 
stirred flame (SF) regime. In the LEM study by IWoosley et al.1 (120091 ). it was found that 
the size of a region, in which a detonation can be triggered, is about 10 km. Since the values 
of £k/rt listed Tabled] are comparable to this size, we assume that Da ~ 1 corresponds to a 
burning zone of size ~ £k/rt- Then T = £k/rt/^(^k/rt) is the turn-over time of the largest 
eddies in the nearly isotropic regime. Since Kolmogorov scaling applies on length scales 
£ < ^k/rt, it follows that T = ^k/rt/ 6- ^ or ^ a > ^' there are broadened flame structures of 
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size smaller than £k/kt- 

In the framework of intermittency theory, the Damkohler number in a region of size 
£ k /rt is given by 



Da 



.2/3 
'K/RT 



,1/3 

'^K/RT ' 



(9) 



Tt, 



where Q K/RT is interpreted as the random dissipation rate on the length scale £k/rt- Hence, 
the requirement Da > 1 for a DDT in the SF regime corresponds to an upper bound on 
Ck/rt: 



^k/rt < £ wsr ■'- 



()2 

'K/RT 



(10) 



If Da becomes less than unity (q k/rt > cwsr), the well stirred reactor (WSR) regime is 
entered. In this regime, the density becomes so low that a detonation is very difficult and, 
eventually, the flames will be quenched. 

On the other hand, if Da is yet too high, the broadened flames produced at the onset of 
distrib uted burning will be too sparse to coalesce into a mixed flame structure extending over 
£k/rt (jWoosley et al.l 120091 ) . For this reason, there is a lower bound on €k/rt corresponding 
to a critical Damkohler number Da cr it: 



6 ^k/r.t > C crit 



ft 

'K/RT 



^ a crit 'nuc 



If a DDT was only constrained by the range of Da, in principle, very small velocity fluc- 
tuations could trigger a delayed detonation at sufficiently low density because of the rapid 
increase of r nuc as the bulk expansion causes the density to decrease. The mechanism of a 
DDT, however, is likely to require velocity fluctu ations that reach a fraction ~ 0.1 of the 
speed of sound (ILisewski et al.ll2000l ; IWoosleyll2007l ). A typical figu re for the minim al velocity 
fluctuation v' min on length scales £ ~ £k/rt ~ 10 km is 500 km s -1 (jWoosleyl 120091 ) . 



Combining the constraints (|TT|) and (|T0|) with the requirement €k/rt > 
the local probability of a DDT becomes 



'K/RT, 



P KWSR > Q K /rt > 

tcrit 



=-erfc 
2 



ln(max(e cri t,6 min )/6) <H K/RT 

2V2 



V2o- £k 



/RT 



-erfc 



ln(max(e WSR ,e min )/e) a t 



V2a iK 



+ 



/RT 



K/RT 



2^2 



(12) 



where o~i K 



/RT 



/j,\n (L/£ k /rt)- For the asymptotic pdf (jSJ) proposed by iRopkel (120071 ). on 
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the other hand, it follows that 



P (ewsR > Q k /rt 



> e 



crit 



exp(a 2 ) 



ai( 



r 



ai 



-a (max(e W SR, e min )£) 



l/oi 
ai/3 



«1 



-a (max(e crit , e min )£) 



ai/3 



(13) 



The function T(x,y) is the incomplete gamma function. 



As explained in Section |3~T| the effective DDT probability depending on Da crit is 

POO 

P D DT(Da crit ) = / P(e W SR > Qk/rt > e crit) pdf(p|G = 0) dp. (14) 

The value of Da cr it has t o be determine d from microphysical studies. A reasonable range 
is 10 < Da crit < 100 (jWoosleyl 120091 ) . Assuming t hat a DDT occurs at a transition 
density near 1 x 10 7 gcm -3 in a region of size 10 km (jWoosley et al.ll2009l ). e cr i t ~ 2.3 x 
10 19 cm 3 s _2 /Da 3 crit . This is, even for Da cr i t ~ 100, several orders of magnitude greater than 
ej ~ 1.1 x 10 10 cm 3 s~ 2 . Also note that e m i n ~ 10 17 cm 3 s~ 2 is higher than e cr i t for Da cr i t = 10. 
For this reason, the probability given by equation (|14p is substantially more constrained 
than the probability Pdbt (Ka) defined by (jHJ). Detailed calculations are presented in the 
following Section. 



4. Number of deflagration-to-detonation transition events 

In the previous Section, we obtained estimates of the effective probability of a DDT 
anywhere close to the flame surface. In order to assess whether a DDT is likely to occur at 
some instant, the expectation value A^ddt of the total number of events has to be calculated. 
Basically, this means that the effective DDT probability has to be multiplied with the number 
of regions of a certain size that can be accommodated within the burning zone. In this regard, 
it is important to account for the actual flame geometry, which becomes extremely folded 
and wrinkled in the course of the deflagration phase. 

Specifying the DDT probability as a function of the Karlovitz number (see formula [8]), 
we face the difficulty that the size £ c of a region, in which a detonation is triggered by 
the Zel'dovich mechanism, strongly varies with the mass density. For this reason, there is 
no simple relation between A/^dt an d PoDT(Ka). In order to calculate A^ddt ; it would be 
necessary to weigh each dp-bin of pdf(p|G = 0) with the number of regions of size £ c (p) 
contained in the volume of burning material at densities from p to p + dp. Since £ c is 
much smaller than the numerical resolution A for p > 10 7 gcm~ 3 , this weighted distribution 
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cannot be inferred from the numerical data. Nevertheless, we can account for the flame 
geometry as far as it is numerically accessible by extrapolating the number of grid cells 
enveloping the flame front, to the number N c of critical regions of size £ c given by the 
mean density of the burning material. If the fractal dimension of the flame front is D, it 
follows that N c ~ ]Va(A/4) d , and, thus, we estimate the total number of DDTs to be 
-ZVddt ~ iV" c PDDT(Ka). Of course, the variation of the mass density over the flame surface 
at a given instant implies that N c might be substantially different from N&(A/£ C ) D even for 
£ c ~ A. Nevertheless, once N c P-£,dt (Ka) exceeds unity by a great margin, a DDT is likely to 
occur. Since the flame surface in type la supernovae is extremely wrinkled on length scales 
much smaller than the integral scale L, we tentatively set D — 3, i. e., the flame front is 
assumed to be space-filling. 



For the DDT criterion of IWoosley et al.l (120091 ). on the other hand, an estimate of A^ddt 



is rather straightforward, because it is assumed that the size of regions, in which a DDT 
can be triggered, is given by £k/rt ~ 10 km. This length scale is independent of the mass 
density, and it is comparable to the size of the grid cells during the pase of the explosion we are 
interested in. Thus, we have AT D dt ~ ^Yk /rt-Pddt (Da CT it ) ; where N K/RT ~ N A (A/£ K/RT ) D 
and PDDT(D a crit) is given by equation (TH|) . 



5. Numerical results 

We interpolated eg from the numerical values of the laminar flame speed and the flame 



width (defined by the temperature profile) listed in Table 1 of I Woosley et al.l (120091 ) as a 
function of the mass density. Table 3 in the same article specifies values of the nuclear time 
scale T nuc in the WSR regime for various densities. To estimate the c ritical length scale £ c , 



we made use of the compilation of data in Table 1 of |Pan et al.l (120081 ). Since there are only 



few values for different mass densities, the numerical evaluation of £ c is rather uncertain. 

The resulting probability P(e£ c > eg) as function of the mass density for Ka = 1 is shown 
in the Figure [2]for t = 0.6, 0.7, 0.8 and 0.9 seconds. Except for the latest instant of time, the 
two solid curves were obtained by substitution of the parameters /x and e inferred from the 
angular structure functions and the radial structure functions (see Tabled]). The differences 
between both cases are small. This supports our assumption that the log-normal distribution 
can be determined on the basis of the the parameters /i ang and e ang as well as /i ra( j and e ra d- 
This is important for t = 0.9 seconds, where the regime of isotropic turbulence (£ < £k/rt) is 
numerically unresolved, and we have only the parameters of the angular structure functions 
available. We find that P(ee c > eg) ~ 1 for p less than about 2 x 10 7 gcm -3 The probabilities 
assuming the time-independent models A, C and E (corresponding to = 10 16 , 10 14 and 
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10 12 cm 2 s 3 , respectively) of |Pan et al.l (120081 ) are indicated by the dashed lines in Figure [2j 



Our data fall in between models C and E. The steeper drop of the probability toward higher 
mass density in compariso n to the models of Pa n et al. is a consequence of the calculation of 



e fl using the recent data by lWoosley et al.l (120091 ) . Since it is considered to be more likely tha t 



a DDT occurs once the Karlovitz number becomes greater than 10 (jWoosley et al.ll2009l ). 
we also evaluated equation ([8]) for Ka 2 = 10, 100 and 1000. The resulting probabilities 
P(q c > Ka 2 Cfi) are plotted as functions of the mass density in Figure [3] for t = 0.9 seconds. 
As one see, the range of densities for which P(q c > Ka 2 eg) is about unity decreases roughly 
by a factor of 2 as Ka 2 is raised from 1 to 1000. 

For the computation of the effective probability P DDT (Ka) as defined by equation ([7]), 
we substituted the probability density functions of the mass density constrained to the flame 
surface, pdf (p\G = 0), which are plotted in FigurelU Comparing Figures [2] and HI it is evident 
that pdf(p|G = 0) significantly overlaps with P(ei c > eg) as a function of density for t > 0.8 
seconds only. The product P(ee c > efl)pdf(p|G = 0) is shown in Figure [5] for t = 0.8 and 
0.9 seconds, and the effective probabilities obtained by integrating the probability densities 
are listed in Table [3] At time earlier than 0.8 seconds, the resulting probability of a DDT is 
low. For t = 0.8 seconds, Pddt would exceed 50 %, if a DDT was triggered right at the onset 
of distributed burning (Ka=l). We emphasize that this is a highly unrealistic assumption. 
However, effective probabilities higher than 50 % result for all investigated values of Ka at 
time t = 0.9 seconds. 

Because of the huge number of regions of size £ c that can be accommodated within the 
flames, an order-of-magnitude estimate of N c as outlined in Section H] implies that Addt 
exceeds unity by many orders of magnitude already at t = 0.6 seconds. At this time, the 
mass density at the flame front is of the order 10 8 gcm~ 3 (see Fig. H]), and £ c ~ 10~ 3 km for 
this density. The number of grid cells of size A « 4km is ~ 10 6 . Hence, N c ~ 10 17 for 
D = 3, and with P DDT (Ka = 1) ~ 10~ 7 (see Table |3J, it follows that A^ DDT ~ 10 10 . For 
Ka = 10, PoDT(Ka) decreases by a few orders of magnitude, but A^ddt is still much greater 
than unity. Even for D = 2, N c ~ 10 13 and Addt ^ 1 for Ka < 10. Thus, if the conditions 
for a DDT would be solely constrained by the Karlovitz number, delayed detonations could 
easily occur at early stages of the explosion. 

The graphs of the l ocal probabilit y -Pfew sR > Q k /rt > e crit) according to the DDT 



constraints proposed by IWoosley et al.l (120091 ) are shown in Figure [6] for t = 0.8 and 0.9 
seconds. For comparison, P(ee c > Ka 2 6fl), where Ka = 10 3 / 2 , is also shown. As Da cri t 
decreases, the range of densities for which P(ewsR > Q k /rt > e c"0 ~ 1 becomes increasingly 
narrow. In contrast to P(ee c > Ka 2 6fi), there is also a cutoff toward lower densities, which 
accounts for the fact that burning and, consequently, a DDT cannot occur at arbitrarily 
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low densities. However, the shape of this cutoff is largely qualitative, because it is based 
on extrapolations of the microphysical parameters in the range of densities lower than 0.6 x 
10 7 gcm -3 . 

Calculating the DDT probability from equation (fl4j) for critical Damkohler numbers 
in the ranee from 10 1/3 to 100, results in very small probabilities for t < 0.7 seconds (see 
Table H]). From iV DDT ~ A^/RT-PoDTlDacrit), it follows that the expectation value of iV DDT 
is less than unity, except for Da crit > 100. For t = 0.8 and 0.9 seconds, on the other 
hand, the values of PDDT(Dacrit) are much higher and A^ddt 3> 1 for all critical Damkohler 
numbers if we assume that there is no bound on the minimal velocity fluctuation (i. e., 
e m in = 0). Consequently, a delayed detonation would occur almost certainly between 0.7 and 
0.8 seconds after ignition. The constrained probability density functions -P(ewsR > Q K /rt > 
e C rit) pdf (p\G = 0) plotted in Figure [7] allow us to assess the possible range of transition 
densities subject to the condition A^ddt > 1- Depending on the value of Da crit , transition 
densities ranging from 0.5 x 10 7 to 1.2 x 10 7 gem for t = 0.8 seconds are most likely. 
For Da cr it = 10 a transition density around 0.7 x 10 7 gcm -3 would be preferred. Since 
jVddt ~ 10 4 , the transition density might be closer to 10 7 gcm -3 though. It is reassuring 
to note that our DDT densities are comparable to the ones determined from fits of ID 
delayed detonation models to observational data. However, an exact match is not expected 
and would not necessarily lead to agreement of multi-dimensional models with observations 
because the detonation wave does not propagate spherically outwards. 

However, the probability of a DDT drops dramatically if the minimal velocity to trigger 
a detonation, v' min , is greater than about 100 kms -1 . In Table [5j the results for v' min = 
0, 200 and 500 kms -1 are compared for different values of the intermittency parameter [i 
at time t = 0.9 seconds. As on e can see , dela yed detonations are virtually excluded if 



^min = 500 kms 1 , as proposed by IWoosleyl (120071 ) . Assuming a smaller minimal velocity of 



200 kms -1 , the results turn out to be extremely sensitive on the intermittency parameter 
/i, while only small variation is found for v' min = Okms -1 . The reason becomes apparent 
from the dependence of the log-normal distribution on the intermittency parameter. In 
Figure [H log-normal probability density functions [1] for £ = 10 km, fi = 0.05, . . . , 0.2, and 
the parameters L and e determined by the angular structure functions at time t = 0.9 seconds 
are plotted. Assuming Da crit = 10, the thresholds e crit defined by equation ffTTl) are indicated 
for several different mass densities, which correspond to cumulative mass fractions of 20, 30, 
and 40 % of burning material. If we solely consider the constraint ewsR > Q > fcrit, the 
intermittency parameter has only a small influence, because the cumulative probability is 
mainly determined by the central parts of the distributions. In the case that q must further 
satisfy, say, q > e m i n = 8 x 10 15 cm 3 s~ 2 corresponding to v' min = 200 kms -1 , the cumulative 
probabilities are given by the right wings of the log-normal pdfs. Depending on the value 
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of \x 1 there are substantial variations of the effective probabilities, as different portions of 
the wings contribute. If e m j n exceeds about 10 17 cm 3 s -2 , the probability of a DDT becomes 
virtually zero. 



The asymptotic tails of the Q-pdfs calculated by iRopkd (120071 ). on the other hand, 
permit delayed detonations even for v' min = 500 kms -1 . Table [6] summerizes the results for 
-Pddt and -/V DDT at t = 0.7, 0.8 and 0.9 seconds. For v' min = 500 kms -1 , DDTs are likely to 
occur between 0.7 and 0.8 seconds. The constrained probability density function at t = 0.8 
seconds is plotted in the left panel of Fig. [9J The transition density for Da crit = 10 is 
close to 1.0 x 10 7 gcm -3 . After 0.8 seconds the probabilities begin to decline because of the 
decreasing turbulence intensity (see Tableland the right panel of Fig. [9]). This behavior is 
in marked contrast to the predictions based on lognormal distributions, for which the DDT 
probabilities at t — 0.9 seconds are higher than for 0.8 seconds (see Table H] and Fig. [7]). 
Therefore, there appears to be a narrow time window, in which DDTs are possible. In Fig.[T0| 
plots of pdf DDT (Da crit ) for t = 0.8 seconds are shown for v' min = 200 and 1000 kms -1 . It is 
obvious that the value of v ' min has a huge impact on the probability density functions. As 
a consequence, delayed detonations could be set off earlier than 0.7 seconds after ignition 
if ^min was about 200 kms -1 or less, whereas the DDT probabilities would be marginal for 
v' miQ = 1000 kms -1 . 



6. Conclusion 



We investigated the intermit tency propert i es of tu rbulence in the n umerical simulation 
of a thermonuclear supernova by iRopke et al.l (120071 ). iPan et al.1 (120081 ) proposed that the 
probability of entering the distributed burning regime can be computed from the log-normal 
model for intermittent turbulence, and from this probability the incidence of a DDT can 
be inferred. Evaluating the characteristic scales of turbulence and the probability density 
functions of the mass density in the vicinity of the flame front at different instants, we 
calculated the probability of a DDT for various Karlovitz numbers. We also i nvestigated the 



i nfluen ce of more restrictive criteria following from the numerical studies by IWoosley et al. 
(l2009h . 



Assuming that a detonation can be triggered after the onset of distributed burning, our 
calculations indicate that a delayed deton ation would commence at an early stage of the 
explosion. In contrast to IPan et al.l (120081 ). where a spherical flame is assumed, the highly 
wrinkled and folded flame front in the numerical simulation greatly increases the number of 
regions of critical size. Although we cannot precisely determine the to tal number of regions 
of size £ c , because the critical size greatly varies with the mass density (INiemeyer fc Woosley 
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19971 ). and £ c <C A, where A is the grid resolution, it appears that TVddt greatly exceeds 
unity already 0.6 seconds after ignition. Then the transitio n density would be signifcantly 
higher than 10 7 gcm -3 . This possibility was pointed out by lPan et all (120081 ). 



If DDT c onditions are constra ined by an interval of Damkohler numbers in the stirred 
flame regime (jWoosley et al.ll2009i ). a DDT earlier than 0.7 seconds after ignition can be 
excluded. Assuming a typical value of the Damkohler number, for which a DDT can be 
triggered, a delayed detonation could occur around t = 0.8 seconds, and the transition 
density would be close to 10 7 gcm -3 . However, assuming a log-normal distribution, the 
typcial velocity fluctuations would b e implausibly low. If turbulent velocity fluctuations 
greater than 500 kms -1 are required (jWoosleyl |2009i ) . then non- log- normal distributions are 
vital for delayed detonations. 

One of the problems of determining log-normal distributions from scaling exponents 
of structure functions is that it is difficult to calculate higher-order two-point statistics. 
Consequently, the computed scaling exponents are sufficient to constrain the central part of 
the distribution only, whereas the tails are based on an extrapolation. A further caveat is that 
we are not able to separate ash from fuel regions, because two-point statistics of turbulence 
can only be computed in convex regions . For D DTs, however, turbulence in fuel close to the 
flame front is signicant. Indeed, iRopkd (120071 ) showed that the tail of the pdf of turbulent 
velocity fluctuations on a length scale of about 10 km in the vincity of the flame front is 
much flatter than what is expected on the basis of log-normal models. With this distribution, 
DDTs between 0.7 and 0.8 seconds after ignition at a transition density close to 10 7 g cm -3 are 
definitely possible. Of course, the detonation time also depe nds on the ignition scenario and 
the subsequent evolution of the deflagration (see, for example, ISchmidt fc Niemeveri[2006l l In 
any case, it is crucial to settle the question of the relevant distribution of turbulent velocity 
fluctuations in the future. Given a certain distribution, the minimal velocity fluctuation 
that is necessary to trigger a DDT is an extremely important parameter. Depending on this 
threshold, delayed detonations might be theoretically confirmed or excluded. Consequently, 
further studies of the microphysics, particularly in the range of densities below 10 7 gcm -3 , 
will be essential for a more accurate evaluation of the DDT probability in type la supernovae. 

Based on the insights provided by such calculations, it may be possible to devise an 
algorithm that determines the local probability of a DDT in large-scale simulations of ther- 
monunclear supernovae. Since the processes causing a DDT occur on scales t hat are difficult 



to res olve in such simulations, a subgrid scale model is likely to play some role (ISchmidt et al. 



20061 ). Apart from this, the p ropagation of t he burning zone has to be treated beyond the 
onset of distributed burning. ISchmidtl (120071 ) proposed that the level set te chnique can be 
extended at least into the broken-reaction-zones regime (IKim &: Menonll2000l ). provided that 
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the burning time scale does not exceed some fraction of the eddy turn-over time correspond- 
ing to the numerically unresolved velocity fluctuations (i. e., the ratio of the grid cell size to 
the square root of the specific subgrid scale turbulence energy). This approach also calls for 
further microphysical studies. An entirely different approach might be the use of adaptive 
mesh refinement an d the in s i tu cal culation of the processes triggering a DDT, for instance, 
by means of LEM (IKersteinl 1 19911 ). Once all numerical challenges are met, quantitative 
theoretical arguments in favor or against delayed detonations as an explanation for type la 
supernovae will be within reach. 
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Fig. 1. — Scaling exponents of the angular (left) and radial (right) structure functions of 
order p < 6 at time t = 0.7 s. For comparison, the predictions of the Kolmogorov theory 
(dot-dashed line), the log-normal model (blue dashed line in the online version) with p = 0.2 
and the She-Leveque model (red dashed line in the online version) are plotted. The solid 
lines show the log-normal and log-Poisson fit functions, which are nearly coinciding. For the 
log-Poisson model, we set the codimension C = 2. The corresponding fit parameters are 
listed in Table CD 



-22 - 




Fig. 2. — Dependence of the probability that the rate of dissipation q c in a region of critical 
size £ c exceeds the threshold on the mass density for different stages of the explosion. For 
each instant, the results inferred from angular and radial struc ture functions (so lid lines) are 
plotted together with three different models (A, C and E) of |Pan et al.l (120081 ). The mass 
density is specified in units of 10 9 gcm~ 3 . 
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Fig. 3. — Variation of P [e^ c > Ka 2 efl) with the Karlovitz number Ka for t = 0.9 seconds. 
From right to left, the curves correspond to Ka 2 = 1, 10, 100 and 1000. The mass density is 
specified in units of 10 9 gcm -3 . 
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Fig. 4. — Probability density functions of the mass density constrained to the flame front at 
different instants. The mass density is specified in units of 10 9 gcm -3 . 
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Fig. 5. — Constrained probability density functions for triggering a DDT by velocity fluctua- 
tions satisfying q c > Ka 2 ea for t = 0.8 (left) and 0.9 (right) seconds. For Ka 2 increasing from 
1 to 1000 by factors of 10, the distribution becomes increasingly narrow. The calculation is 
based on the parameters \i and e obtained from the angular structure functions. The mass 
density is specified in units of 10 9 gcm -3 . 
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Fig. 6. — Variation of P (ewsR > Q > e crit) with the critical Damkohler number, Da crit , for 
I = 10 km at time t = 0.9 seconds. The distribution becomes increasingly narrow as Da^. rit 
decreases from 10 6 to 10 in order-of-magnitude steps. For comparison, P (ee c > 1000efl) is 
shown as dashed curve. The mass density is specified in units of 10 9 gcm -3 . 
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Fig. 7. — Constrained probability density functions for triggering a DDT by velocity fluctu- 
ations satisfying ewsR > Q > ^crit for t = 10 km at time t = 0.8 (left) and 0.9 (right) seconds, 
where the solid curves correspond to Da^. rit decreasing from 10 6 to 10 in order-of-magnitude 
steps from right to left. The constrained PDF resulting from the criterion q c > lOOOea is 
indicated by the dashed line. The calculation is based on the parameters \i and e obtained 
from the angular structure functions. The mass density is specified in units of 10 9 gcm -3 . 
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Fig. 8. — Log-normal probability density functions (CQ) of €i in units of 10 15 cm 3 s~ 2 for 
I = 10 km at time t = 0.9 seconds. The mean rate of dissipation e and the integral scale 
L are determined by the third-order angular structure functions (see Table [1]). The four 
plotted pdfs follow for mu = 0.05 (most narrow distribution), 0.1, 0.15, and 0.2 (widest 
distribution). The tails beyond 5<T£ are dashed. The vertical lines indicate the values of e crit 
defined by equation fllip for different mass densities, assuming that the critical Damkohler 
number is equal to 10. The percentages indicate the volume fractions of matter at densities 
less than the mass densities chosen for the evaluation of e cr i t . 
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Fig. 9. — The same probability density functions as in Figure 0, however, with the local 
DDT probability given by equation ( fl3l) in place of (1121) . and v' min = 500 kms" 1 . Note that 
the ordinate scale is different from Figure [71 
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Fig. 10. — Variation of the probability density functions shown in the left panel of Fig. [9] 
for v' min = 200 kms -1 (left) and v' min = 1000 kms -1 (right). Note that the different ordinate 
scales. 
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t[s] 


L 


Ik/rt 


fang 


^rad 


/^ang 


/^rad 


/^ang 


A-ad 


A ang 




0.5 


516 


15.4 


0.0332 


0.0634 


0.096 


0.139 


0.783 


0.738 


0.434 


0.523 


0.6 


799 


14.4 


0.0819 


0.0972 


0.128 


0.134 


0.750 


0.744 


0.501 


0.512 


0.7 


1023 


14.3 


0.166 


0.177 


0.111 


0.121 


0.766 


0.755 


0.468 


0.489 


0.8 


1578 


12.9 


0.335 


0.274 


0.140 


0.128 


0.738 


0.749 


0.525 


0.503 


0.9 


2324 




0.551 




0.132 




0.745 




0.510 





Table 1: Parameters of log-normal and log-Poisson model fits to the angular and radial 
scaling exponents at different instants. The radial scalings apply to the subrange £ < £k/rt- 
Also specified are the integral scale and the mean rate of dissipation inferred from the radial 
and angular structure functions. 



t [s 


A'' ang 


2 C6,ang 


Std(C6,ang) 


A'rad 


2 C6,rad 


Std(C 6 ,rad) 


0.5 


0.096 


0.099 


0.037 


0.139 


0.144 


0.047 


0.6 


0.128 


0.125 


0.033 


0.134 


0.141 


0.055 


0.7 


0.111 


0.109 


0.028 


0.121 


0.129 


0.057 


0.8 


0.140 


0.141 


0.044 


0.128 


0.134 


0.054 


0.9 


0.132 


0.132 


0.043 









Table 2: Fitted intermittency parameter of the log- normal model and values following from 
the sixth-order scaling exponents of the angular and radial structure functions. Also listed 
are the standard errors of the scaling exponents following from power-law fits to the structure 
functions (see paper I). 



t[s] 


Ka 

p( an s) 

-'DDT 


= 1 

p(rad) 
-'DDT 


Ka 

p( an s) 

-'DDT 


= 10 

p(rad) 
-"DDT 


Ka = 

p( an g) 
-'DDT 


31.6 

p(rad) 
-'DDT 


0.6 


3.0 X 10~ 7 


4.8 X 10~ 7 


~ 


~ 


~ 


~ 


0.7 


0.045 


0.050 


0.0024 


0.0026 


4.2 x 10- 4 


4.8 X 10~ 4 


0.8 


0.52 


0.53 


0.22 


0.22 


0.12 


0.13 


0.9 


0.93 




0.72 




0.59 





Table 3: Effective probability of a DDT inferred from equation (J7J) for different values of the 
Karlovitz number. 
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Da cr j t 


p(ang) 
DDT 


7\r( an g) 
DDT 


p(rad) 
-'DDT 


7\r(rad) 
DDT 


t = 0.7 seconds, N K/RT « 6.1 • 10 5 


2.15 


~ 


~ 


~ 


~ 


4 64 


~ 


rsj 


~ 


~ 


10.0 


~ 


~ 


~ 


~ 


21.5 


1.1 x 10" 


8 7.0 x 10- 3 


rsj 


~ 


46.4 


4.4 x 10" 


-e 2.7 


4.1 x 10" 6 


2.5 


100 


7.6 x 10" 


" 5 4.7 xlO 1 


8.6 x 10~ 5 


5.3 x 10 1 




t = o.t 


I seconds, N K/RT « 1.7 • 10 6 




2.15 


0.0021 


3.9 x 10 3 


0.0022 


3.7 x 10 3 


4.64 


0.0061 


1.1 x 10 4 


0.0061 


1.1 x 10 4 


10.0 


0.013 


2.4 x 10 4 


0.014 


2.4 x 10 4 


21.5 


0.026 


4.8 x 10 4 


0.027 


4.7 x 10 4 


46.4 


0.046 


8.4 x 10 4 


0.047 


8.2 x 10 4 


100 


0.076 


1.4 x 10 5 


0.077 


1.3 x 10 5 


t = 0.9 seconds, N K/RT « 4.1 • 10 6 


2.15 


0.04 


1.8 x 10 5 






4.64 


0.10 


3.9 x 10 5 






10.0 


0.15 


6.3 x 10 5 






21.5 


0.22 


8.8 x 10 5 






46.4 


0.28 


1.2 x 10 6 






100 


0.36 


1.5 x 10 6 







Table 4: Dependence of the effective probability of a DDT and the expectation value of the 
number of DDTs on the critical Damkohler number for several instants. 





v' ■ 

mm 


= km s 1 


v' ■ = 

mm 


200 kms- 1 


v' ■ = 


500 kms" 1 




-Pddt 




-Pddt 


A^DDT 


-Pddt 


iVDDT 


0.05 


0.155 


6.34 x 10 5 


6.4 x 10~ 


9 2.6 x 10~ 2 


~ 


~ 


0.1 


0.154 


6.30 x 10 5 


5.8 x 10- 


6 2.4 x 10 1 


~ 


~ 


0.15 


0.153 


6.25 x 10 5 


5.8 x 10- 


" 5 2.4 x 10 2 


~ 


~ 


0.2 


0.152 


6.20 x 10 5 


1.8 x 10- 


4 7.5 x 10 2 


1.1 x 10- 


9 4.4 x 10- 3 



Table 5: Dependence of the effective probability of a DDT and the expectation value of the 
number of DDTs on the intermittency parameter \x of the log-normal model and the minimal 
velocity fluctuation t> r ' nin for t = 0.9 seconds. 
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v' ■ = 

mm 


200 km s" 1 


V 'mm = 500 km s 1 


v' ■ = 

mm 


1000 km s^ 1 




Da cr jt 


-Pddt 




A^DDT 


-Pddt 


-/Vddt 


-Pddt 




-/Vddt 




t = 0.7 seconds, iV" K /RT « 6.1 ■ 10 5 


2.15 


~ 




~ 


~ 


~ 


~ 




~ 




4.64 


3.7 x 10" 


-9 


2.3 x lO" 3 


2.3 x 10~ 9 


1.4 x 10" 3 


~ 




~ 




10.0 


3.9 x 10" 


-7 


2.4 x 10" 1 


4.2 x 10~ 8 


2.6 x 10- 2 


~ 




~ 




21.5 


5.9 x 10- 


-6 


3.6 


2.1 x 10~ 7 


1.3 x 10" 1 


~ 




~ 




46.4 


3.0 x 10' 


-5 


1.9 x 10 1 


6.6 x 10~ 7 


4.0 x 10" 1 


1.7 x 10- 


-9 


1.0 x 10- 


-3 


100 


9.5 x 10- 


-5 


5.8 x 10 1 


1.8 x 10~ 6 


1.1 


4.5 x 10- 


-9 


2.8 x 10- 


-3 








t = 0.8 


seconds, Nk/rt ~ 1-7 • 10 6 










2.15 


3.5 x 10- 


-4 


6.2 x 10 2 


6.2 x 10~ 6 


1.1 x 10 1 


3.9 x 10- 


-8 


6.8 x 10- 


-2 


4.64 


9.9 x 10" 


-4 


1.7 x 10 3 


1.6 x lO" 5 


2.8 x 10 1 


9.8 x 10- 


-8 


1.7 x 10- 


-1 


10.0 


2.0 x 10- 


-3 


3.5 x 10 3 


3.1 x 10~ 5 


5.3 x 10 1 


1.8 x 10- 


-7 


3.2 x 10- 


-1 


21.5 


3.5 x 10- 


-3 


6.1 x 10 3 


5.1 x 10~ 5 


8.9 x 10 1 


3.0 x 10- 


-7 


5.1 x 10- 


-1 


46.4 


5.6 x 10- 


-3 


9.8 x 10 3 


7.9 x 10~ 5 


1.4 x 10 2 


4.5 x 10- 


-7 


7.8 x 10- 


-1 


100 


8.5 x 10- 


-3 


1.5 x 10 4 


1.2 x 10~ 4 


2.0 x 10 2 


6.3 x 10- 


-7 


1.1 





t = 0.9 seconds, N K/RT w 4.1 ■ 10 ( 



2.15 


6.1 x 10" 


-4 


2.5 x 10 3 


1.2 x 10" 


-6 


4.8 


~ 




~ 




4.64 


1.3 x 10" 


-3 


5.3 x 10 3 


2.4 x 10" 


-6 


9.9 


1.2 x 10- 


-9 


5.0 x 10- 


-3 


10.0 


2.0 x 10" 


-3 


8.3 x 10 3 


3.8 x 10" 


-6 


1.6 x 10 1 


1.9 x 10- 


-9 


7.8 x 10- 


-3 


21.5 


2.8 x 10" 


-3 


1.2 x 10 4 


5.2 x 10" 


-6 


2.1 x 10 1 


2.6 x 10' 


-9 


1.1 x 10- 


-2 


46.4 


3.7 x 10- 


-3 


1.5 x 10 4 


6.7 x 10- 


-6 


2.8 x 10 1 


3.3 x 10- 


-9 


1.4 x 10- 


-2 


100 


4.6 x 10" 


-3 


1.9 x 10 4 


8.3 x 10- 


-6 


3.4 x 10 1 


4.0 x 10- 


-9 


1.6 x 10- 


-2 



Table 6: Dependence of the effective probability of a DDT and the expectation value of the 
number of DDTs on the critical Damkohler number co rresponding to Table H] for a non-log- 
normal distribution of turbulent velocity fluctuations (IRopkel 120071 ). 



